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Bound states of the pion-nucleus system are investigated in momentum space using a micro- 
scopic optical potential with inherent energy dependences and nonlocahties arising from elementary 
potential models. The wave equation and computational techniques are tested, shallow levels are 
calculated, and a comparison is made with experiment. Deep levels are found to be deeper and 
broader than those in other investigations. 
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I. INTRODUCTION 

rn ■ 

I A complex potential is sometimes used in quantum mechanics to model a physical system which is coupled to 
another system or channel. This approach is used, for example, in analyzing exotic atoms where an optical potential 
is added to the Coulomb potential to describe the exotic particle's strong interaction and eventual annihilation with 
f^. the nucleus. The normalizable states for such systems are "quasibound" since they decay exponentially in time with 
^ ' a lifetime proportional to the inverse of the imaginary part of the binding energy. 

Because strong channel coupling is non perturbative, its phenomenology is often non intuitive. For example, if 
'- a potential has sufficiently large imaginary (absorptive) part, an increase in the absorptive strength may lead to 
[ a decrease in the level's width [Q. Here the increase in the imaginary potential leads to further absorption of the 
^-H ' wavefunction from the nuclear interior, and the corresponding decrease in the wavefunction's overlap with the nucleus 
J> [ leads to a decrease in the level's width. In these cases, even if the real part of the optical potential is attractive, the 
1 large imaginary part acts as a repulsive potential by excluding the wavefunction from the nucleus. 

Recently, theoretical studies by Friedman and Scoff Q, Hirenzaki et al. and Toki and Yamazaki [QH] have 
pointed out that for deeply bound levels in heavy nuclei, the pion optical potential might be sufficiently absorptive 
' for there to exist states of the type just described. Furthermore, these quasibound states would also be "quasistable" 
I (and thus observable) since their widths would be smaller than the surrounding energy-level spacings (F ~ AE/10). 
^\ • For example, ^"^Pb would have Is and 2p levels bound by ^ 7 and ~ 5 MeV and widths ~ 0.5 MeV — in contrast to 
] the ~ 20 MeV widths expected if lowest-order perturbation theory were valid. Furthermore, although the pion is close 
> . enough for it to get bound up with the nucleus, the absorption so excludes the wavefunction that the states become 
" hybrid Coulomb-nuclear ones with the dominant binding arising from the Coulomb potential, and the probability 
O density peaking just outside the nuclear surface. 

Observing deeply bound pion states in heavy nuclei is difficult because the X-ray cascade ends with the pion being 
absorbed before it reaches the deep levels. Transitions to the Is level are seen experimentally only for Z < 12 and to 
the 2p level only for Z < 30, and so alternative observational methods have been proposed. A recent experimental 
search using the {n,p) transfer reaction did not record the distinct peaks expected for a Is or 2p transition in Pb 
^ 1^. Apparently, the cross sections are smaller, the widths larger, or the binding energies less certain than predicted. 
d Other suggested experiments include {n,d) pickup 0], (7r~,7) radiative trapping and {n^,p) pickup 

The calculations of deeply bound pion states generally use a phenomenological, coordinate space, optical potential 
of the Ericson-Ericson form: 

— Att 

VP'ir) — [bop{r) + coV ■ p{r)V + Sop'(r) + CqV • p'{r)V] (1) 

where for clarity we assume here an isoscalar nucleus and leave off the Lorenz-Lorentz factor. The parameters 
6o J Co, -Bo I Slid Co generally are determined by ffts to less-bound energy levels in pionic atoms. Because the pion has 
light mass, relativistic effects are expected to be important, and because the pion has zero spin, the correct dynamical 
equation is assumed to be the Klein-Gordon equation: 

[S-y™-'(r)]'^(r) = [-V2 + {/i + y°P*(r)}'] ^(r) (2) 
[£:_yc°«'(j.)]2^(r) ^ [-V2 + /i2 + 2Aiy°P*(r)] V'W (3) 

In (^ the optical potential is treated as a scalar and added to the mass p and the Coulomb potential is treated as 
the time component of a four-vector and subtracted from the energy E. Because the square of the optical potential 
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introduces technical difficulties, the approximation (||) with no (1/°^*)^ is usually used |^;| 10| (Kwon and Tabakin 
PI have indicated how the quadratic term can be handled in momentum space). 

The investigation we report upon here was stimulated by the question of how much the properties of deeply bound 
pionic states would change if a different, possibly more microscopic, form for the optical potential and wave equation 
were used. In particular, we wondered about using LPOTT, the potential we had previously developed for scattering 
1^,^, but never used for bound states. LPOTT is a theoretical optical potential in momentum space having 
nonlocality, off-shell unitarity, and energy dependence (in addition to its momentum dependence) arising from that 
of a pi-nucleon T matrix derived from a separable potential. Since states bound by ^ 10 MeV are quite far down in 
energy from those at which the constants bo,co,BQ, and Cq were fit, we particularly wondered if LPOTT s inherent 
energy dependence might make a difference. Other factors of interest and difference include finite ttN interaction 
range effects in both the first and second order potentials, covariant kinematics and transformation of scattering 
amplitudes ("angle transformation"), and the use of the relativistic Schrodinger equation. (Since the relativistic 
Schrodinger equation is also used in the definition of the ttN T matrix and in the derivation of the optical potential, 
our framework is consistent, which unfortunately, use of (|l|) in (^) is not.) 

In §^ we describe the equations and computat ional method used. In § 111 A we examine relativistic effects and the 
numerical precision of our c ompu tations, and in § 111 B| make comparisons to the experimental ene rgies a nd widths of 
shallow states. Finally, in § 111 C| we examine the deep states supported by our potential, and in § IIID deduce what 
changes in the potential would be needed for there to exist pure nuclear pion bound states. This in turn leads to 
further understanding of the low energy optical potential. 



II. FORMALISM AND COMPUTATIONAL METHOD 

The momentum space relativistic Schrodinger equation is 

H\^j) - E\ij) (4) 
K{k)Mk) + - Vi{k,k';E)Mk')k' dk' =EMk) (5) 

Jo 

y^ml + k'^ + yjm\ + k'^ - - TTiA = K{k) (6) 

Here K{k) is the relativistic kinetic energy and Vi{k,k ; E) is the energy-dependent pion-nucleus potential after 
partial wave projection. There is no difficulty in momentum space with the square roots in K{k) since they contain 
just numbers not operators. Although the relativistic Schrodinger equation (^-(^ obviously contains relativistic 
kinematics, it is not covariant. However, it is the direct part of the Blankenbecler-Sugar equation, which is the 3D 
reduction of the Bethe-Salpeter equation, a covariant, two-particle equation ||l^ , |l5|] . The Klein-Gordon equation, in 
turn, is relativistic and covariant, but it is only a one-particle equation, and when used as in (^ does not treat the 
optical potential covariantly. 

To handle the most general case of coupled and open channels, we transform (^) to an integral equation involving 
the Green's function (where the appropriate boundary conditions are incorporated) ]T^: 

l^P) =GEVE\ij) (7) 

2 1 r°° 

^'(^) ^ - F T^^i.^ / Vi{k,k';E)Mk')k"dk' (8) 
TT E - K[k) Jq 

Note that for bound states there is no incident wave in (0) and no +ie in Ge- If we consider Ve,K, and G^; as 
operators or their matrix representation, then we see that solution of the eigenvalue problem (^ is equivalent to 
demanding nontrivial solutions of (0), namely, 

(1 - GbFb) = det\l-GEVE\=0 (9) 

We search for solutions of (^) in complex energy space after removing the Coulomb singularity with the Lande' 
subtraction technique. Solving (^ is equivalent to determining the energies of the poles of the 7r-nucleus T matrix 
In contrast, the momentum space formulations of Kwon and Tabakin ||ll[ and Cieply et al. [|l^ directly solve (|^) 
as an eigenvalue problem (Kwon and Tabakin use the Lande' subtraction and Cieply et al. use the Vincent -Phatak 
cutoff). In further contrast, Kalbermann et al. [T^ sidetrack momentum space entirely and include finite range effects 
(but not the full energy dependence of a separable nN T matrix) , by solving an integro-differential wave equation in 
coordinate space. 
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In our calculations, the potential Ve is taken as the sum of electromagnetic plus optical potentials: 



{k\VE\k) = y-"'(<7) + + V^''{k,k) (10) 

y-'i't)--^^ (11) 



,VP,^_ Za^Ppiq) r (2^2 + l)(i2_ 1)1/2 



Here the momentum transfer q = k' — k, y^""' is the Coulomb potential for a nucleus of finite size, Pp{q) is the proton 
form factor, V^^ is the order-Za^ vacuum polarization (Uehling) potential and Ae is the electron Compton 

wavelength. The optical potential V'^^ is evaluated in the impulse and factored approximations, and contains terms 
of first and second order in p: 

vE^'ik , k) = [k , k) + U'i^ [k , k) (13) 
U^,^\k',k) ^ Z{k\po - <r|CL I ^'Po) Ppiq) + N{k',po - g |CL | k,Po) Pn{q) (14) 

= ^Z^ n \ -Pi(cos6>^'jA(g) (15) 
D. {co)^ — --r ^ /y^} , , (16) 

(^'^^^^ 2^.427r)3 [^°(^) gg(.o) ^^"^''^ .g,^(.o) " -"^^ ^'^^ ^'^^ 

^^('^)-^^' (18) 

Here labels the TriV eigenchannel (« = proton, neutron, or isospin, I = orbital, and j = total angular momentum), 
(Jiij is the sign of the potential in that eigenchannel, and k and k' are the initial and final ttTV COM momenta. 
The nuclear form factors Pp^niq) in ( p^ ) and ( p^ are Fourier transforms of Woods-Saxon densities |12 . To aide 



comparison with other optical potentials, we take the size parameters as the standard set used by Friedman et al. |2C] 
with extension to Pb by use of the charge parameters from De Vries et al. |2l) . The parameters are given in Tablc| . 

For simplicity of calculation we have used the factored form of C/'-^-' (p^, which clearly does not average over the 
nucleons' Fermi momenta. We can remove this limitation, yet the computation would be longer, more complicated, 
and possibly less reliable as we extrapolate deep down in energy. We do include nucleon recoil in part, by making the 
optimal choice for the momentum of the struck nucleon [l^ , 

k (A-l)q 

This po is optimal in producing the best factored approximation and the most physical off-energy-shell kinematics. 
We include further recoil and Fermi motion effects in our use of the three-body subenergy losb as the TriV energy in 
( |l^ ) . This choice arises from a three-body model for the optical potential [|3| in which there is a pion of momentum 
fc, a nucleon of momentum po, and a core of momentum P = —{k+p + po). The nN T matrix is then evaluated at 
an energy which is the magnitude of the difference in energies of the irN pair and the core: 

iolBiE) = {k>: + k'X- ^ + rn, + mA- Ea-i{P) - Eef - P^ (20) 

Here Eb is an effective core-nucleon binding energy we take as 20 MeV. Changing Eb to 10 MeV would increase our 
final ri^(Pb) by ~ 20%, and Ere{Ph) by ~ 1%. 

Pauli blocking is included by modifying the t:N T matrix for a nucleon embedded in a Fermi sea with Fermi 
momentum Pp appropriate to the nucleus in question [ p^ . Specifically we modify the D functions in (|l6|) with an 
angle-averaged Pauli operator Qq: 
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Qo, otherwise 



i^p+^r 



p2 



A^Pk 
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(21) 
(22) 
(23) 



The two-body T matrices and momenta in ( [Til ) are related to off-shell ones in the ttN COM by a Lorentz covariant 

prescription IT^ which determines 7, k, and k . As indicated by the guj product in (|lj), we include the finite range of 
the ttTV interaction by using a separable potential which, in turn, was fit to TriV scattering data. Note, the explicit 
energy dependence given by the solution of the ttN T matrix in (|4|)-(^5|) is important for the present study of deeply 
bound states because it permits a physically motivated evaluation of optical potential far below the nN threshold. 
Since our model has finite range ttN interactions, inclusion of a Lorentz-Lorenz term is inappropriate [|l3|,^ . 

The second order potential C/*-^^ includes contributions from processes such as virtual nuclear excitation and pion 
annihilation (true absorption). It is particularly important for pionic atoms since its imaginary part provides ab- 
sorption. The form we use for U^"^^ in (17) is a generalization of the separable form used for U'^^^ (we take kq — > 
for subthreshold extrapolation). Unfortunately, while tremendous theoretical and experimental progress has been 
made in understanding the contribution of true pion absorption to the pi-nucleus interaction, that progress is not 
far enough along to permit a determination of C/^^^ with enough accuracy for present purposes. Consequently, the 
quantities Bp and Cq are fit to the Is and 2p levels' shifts and widths in pionic ^^O and '^"Ca. In addition, while for 
scattering we modeled the energy dependence of and Co as that of the tt+Z? — > pp cross section, we have not been 
able to determine a reliable, analytic continuation for subthreshold energies using the energy dependence of either 
a{'K^D pp) or the parameters of Cieply et al. |l^. While our estimates indicate that the subthreshold energy 
dependence is weak, our ignorance of it weakens our model. 

Since our solution of (^) determines the complex eigenenergies but not eigenfunctions, once the eigenenergies are 
known we solve for the momentum space wavefunctions by using inverse iteration [Tl| ] on the Schrodinger equation 



The wavefunction normalization is then 

1 = 

1 = 

= 



iPfik) k^dk 

[(Re Vj(fc))2-(Im Vj(fc))2] k^dk 



Re Tpiik)!!!! ipi{k)k'^dk 



(24) 
(25) 

(26) 
(27) 
(28) 



Note that the normalization condition (|2^) , which follows from the Gamow state formalism of Hernandez and Mon- 
dragon |p2| , normalizes the square of the complex wave function rather than the square of the modulus. This in turn 
requires the overlap of Re ipiik) and Im ipi{k) to vanish (|2^), which often introduces an oscillation in Re '4>i{r). Once 
the momentum-space wavefunction is known, it is transformed to coordinate space via: 



_ ui{r) 



r I - 

TT 



1/2 



Mk)jiikr) k^dk, 



(29) 



III. RESULTS AND DISCUSSIONS 



A. Relativistic Effects and Numerical Precision 



Relativistic effects should be im 
relativistic Schrodinger equation ( 



jortant for deeply bound states in heavy nuclei. While our dynamical equation is the 
|) , in contrast to the more conventional approximate Klein-Gordon equation (||) , it 
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is reasonable to wonder if either equation is adequate for calculating such states. We examine the point Coulomb case 
since then we know the analytic solutions to the Schrodinger equation (the Bohr energies) and to the Klein-Gordon 



equation |15y2 



= , (30) 

{1 + {ZaYI{n - ; - 1 + + 1 )2 _ (^a)2]2}i/2 



■mc'^{Za)'^ mc^{Za)'^ f n 3 



O(a^) (31) 



To compare with the two-particle Schrodinger theory, we use a reduced mass fi in the second term of ( ^l| ) and the pion 
mass in the first term. We examine in Tables || and III both tt — ''"Ca and tt — ^"^Pb because Za is greater than 
i for Pb which means the Klein-Gordon, point-Coulomb solution ( |30| ) is pathological for Z = (the finite-nucleus 
case is not). 

Columns 2 and 3 in Table ^ indicate our numerical precision by comparing the computed and analytic results for 
the nonrelativistic Schrodinger equation (we know of no analytic results for the relativistic Schrodinger equation). We 
see that the numerical results are good out to the fifth decimal place, and that the computed 2s and 2p levels are 
degenerate (as well they should be) within numerical precision — even though the partial-wave potentials are quite 
different. Columns 3 and 4 in Table || and columns 2 and 3 in Table II] indicate that relativistic effects are significant: 



2% for 40Ca-ls^O% for ^ospb-ls, and 48% for 208pb-2s. In contrast, columns 5 and 6 in Table y and columns 3 
and 4 in Table |l| indicate that both the relativistic Schrodinger equation and the Klein-Gordon equation using the 
reduced mass, remove the degeneracy with I and include similar enough relativistic corrections for the differences 
to lie in the fourth significant figure. (The effect would be bigger for s states in Pb — if only they existed for the 
Klein-Gordon equation.) 

Because these comparisons indicate that the inclusion of terms of higher order in the potential are important, 
and because the deeply bound states are hybrid nuclear-Coulomb, we conclude that ignoring the (]/°p*)2 term in 
the Klein-Gordon equation (|^) may have a significant (probably several percent) effects upon the answers. This is 
consistent with the smaller effect found by Kwon and Tabakin |jll| for the lower energy, and less relativistic, 3d — 2p 
transition in '^^S. 



B. Experimental Shifts of Shallow States 

The only adjustable parameters in our theory are the annihilation strengths Bq and Cq. Because our ultimate 
interest is the Is and 2p levels (and because we want to avoid anomalies associated with the 3c? levels [24|), we fit 
the experimental energies and widths for the Is and 2p levels in pionic ^^O |^,^ and Is level in ""^Ca [27,g^, and 
determined 

Bo = (-0.074 + 0.067^)m-^ Co = (0.051 + 0m9i)m~^ (32) 

Since J/'^' provides the absorption in our model, the imaginary parts of Bq and Co should have the most similarity 
to the work of others. The value for Im Co in ( p2|) is similar to those found in the coordinate space models, while the 
value for Bq is larger, although the momentum-space pionic atom calculation of Cieply et al. are more relevant. 
They fit the Is levels of ^^c and ^^O, and the 2p levels of ^"^S and '"'Ca, and find Bo = (-0.093 + iOM2)m^^, Co = 
(—0.125 + 0.090i)m~^. The differences probably reflect the differences in U^^^ and in the data fit. 

We do not expect a theoretical potential to give the same level of agreement as a phenonenological potential 
whose many parameters have been determined in a global search. Nevertheless, in Table |^ and Figure |l| we show a 
comparison between the strong interaction shifts predicted by our potential (the -|-'s) and various pionic atom data 
p5|-p0|. We see very good agreement for the heavier nuclei, good agreement for the lighter nuclei, and an incongruently 
large deviation for (recall, we fit Bo and Co to ^^O and "^^Ca). In general, our level of agreement is comparable 



to the finite range, momentum space calculations of Cieply et al. |17| — even though the potentials and calculational 
framework differ. 
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C. Deep States 

In Table ^ and Figure ^ we show the energies and widths of the deeply bound 7r~^°*Pb states. The subscript EM 
denotes the inclusion of only electromagnetic interactions (Coulomb + vacuum polarization), tot the combined EM 
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plus optical potentials, and rl, r2, and r3 the r-space calculations of Toki et al. g, Nieves et al. [|0[, and of Konijn 
et al. 13^. We see that our calculated binding energies are slightly larger (~ 7%) than those of Toki et al., Nieves 
et al., and Konijn et al. |3l], and our widths are significantly larger: for the s states, ~ 20% relative to Toki et al. 
H and ~ 200% relative to Nieves et al. for the p and d states, 50% relative to Toki et al. and 200 - 300% 
relative to Nieves et al.. 

Because the differences increase with binding, the cause appears to be the energy dependences of the potentials, the 
different wave equations used, and the different extrapolations in the complex energy plane. Since the differences are 
much greater for T than Re E, the different annihilation potential has much to do with the differences (we use a finite 
range model, the others do not). We clearly are finding a much greater optical model dependence than (0.4, 15)% in 
(Re E, r) found by Nieves et al. (possibly their 7 — 75% phenomenological renormalization of strengths reduced the 
model dependence they found). In our calculations, the vacuum polarization increases the binding by 0.5% and the 
width by 2%, which may be significant, but not compared to the model dependence. 

To further understand the physics of the nonlocal optical potential and test the momentum space calculation, we 
have studied the wave functions of the deeply bound states. The momentum space wavefunctions are calculated via 
( psi ) and the coordinate space ones via the Bessel transform (^9|). In the top of Figure || we display the squared 
modulus of the momentum space wavefunction on a logarithmic scale, and in the bottom of the figure we display the 
coordinate space wavefunction ipi{r), both for the Is state in Pb. We see in the comparison to the pure Coulomb 
wavefunction, that the optical potential introduces structures into the momentum space wavefunction similar to those 
present in the potential. These wavefunctions are similar in both shape and magnitude to those of Toki et al. and 
should yield similar predictions when used in a DWIA calculation of the formation rate of the state. 

Examination of the r-space wavefunction in the bottom of Figure |^, shows that the strong interaction repels the 
wavefunction out from the nuclear interior and that there is a node in Re V'(^) even though this is a Is state. This 
node is a consequence of the existence of Im ipi{r) and the complex normalization constraint (p6|)-(p8|). The repulsion 
is also evident in Figure ^ where we plot the probability density |3^ \ui{r)\'^ = \ripi{r)\'^ when there is only the EM 
potential (dashed curves) and when there is the EM plus nuclear potentials (solid curves). (The node is now hidden 
by the multiplication.) We see in this figure the repulsion for the Is, 2p and 2s states in Pb, while for 3d and higher 
states there is an attraction. Similar results were found by Toki et al. p]. 

While some repulsion (and the narrow widths) arise from the strong absorption of the wavefunction near the origin, 
the p wave repulsion occurs only for heavier nuclei. To understand this, we go back to the optical potential (|l^), and 
note that for nuclei which have Z ^ TV, the Zt^^ + NV^ term in the theoretical optical potential is attractive in p 
wave — as confirmed by the p level shifts in Table IV, Yet for nuclei which have N » Z like Pb, the i'^" term 
dominates and the p wave potential becomes repulsive. In the left part of Figure H we present V{p,p), the diagonal 
parts of the momentum space potentials for the Is state in Pb and the 2p state in "Ca. The potential in this figure 
is clearly repulsive for the Is state in Pb and attractive for the 2p state in Ca (although not shown, the 2p potential 
in Pb is repulsive as expected). We also see in both cases that the imaginary part of the potential is absorptive. 



D. Nuclear Bound States 



We have found in our numerical searches that the optical potential is attractive enough to produce broad pion- 
nucleus resonances, but not attractive enough to produce pure nuclear bound states. We wondered how much stronger 
the optical potential is needed to support nuclear states, and so undertook a series of computations in which we 
progressively increased the strength of the optical potential while keeping the Coulomb potential constant. We found 
that the hybrid Coulomb-nuclear states eventually became nuclear ones, explicitly, for the range of strengths 

the complex binding energies found are: 

- (0.4 + O.OOH) MeV <Eca<-{^ + 5i) MeV (34) 
-(7 + QAi) MeV < Epb < -(22 + ISi) MeV (35) 

In Figure ^ we present the probability density |Q |u;(r)p for the Is state in 7r-Pb for varying optical potential 
strengths. We see the hybrid state slowly moving closer to the nucleus, until a transition takes place and at fifteen- 
fold strength the state is localized completely within the nucleus. When the optical potential is this strong, the state 
remain bound, for both Pb and Ca, even if the Coulomb potential is switched off. 



The increase in binding and eventual nuclear state is expected for p states in Ca where the level shifts in Table IV 



are attractive. We were surprised, however, to find it for s states where the level shift is repulsive. At least for local 
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potentials, increasing the strength of the potential does not change its effect from repulsive to attractive. If we go 
back and look in the left part of Figure ^, we see that Im Vop{p,p) always has the sign expected for an absorptive 
potential and that Re Vop{p,p) appears attractive for the 2p state in Ca and repulsive for the Is state in Pb. So this 
does not explain how we can get a Is nuclear state. 

In an attempt to unravel the physics behind our nonlocal potential's Is bound states, we made use of our knowledge 
of the radial wavefunction u;(r) = ripi{r) obtained via ( ^5| ) and ([29|), to deduce an equivalent and local coordinate 
space potential. We assumed there is a Vequiv (r) which when used in a nonrelativistic Schrodinger equation produces 
this same ui{r): 

^equiv (^) 



(36) 



We evaluate the second derivative in ( pq) numerically using the central difference algorithm, and obtain the equivalent 
potential shown in the right hand side of Figure ||. While we do not expect such a Vequiv i^) to be the most physical of 
potentials, it is revealing, and the fact that the imaginary part of Vequiv ('') is always absorptive lends some confidence. 

We see in the lower right-hand corner of Figure |^ that the real part of the Vequiv (r) for the normal strength 2p 
state in Ca has a dominant, long range attraction — but also a short range repulsion. And since ui{r) has no nodes, 
there is no division by zero in (^) to cause a sign change in Vequiv (r). Yet we also find that the character of the 
potential changes dramatically as the nuclear bound state forms, with the potential acquiring an oscillation in the 
surface — much like the one expected from the V • p{r)V term in the Ericson-Ericson potential. We find the bound 
state getting more bound, the wavefunction developing another node, and this showing up in the effective potential. 

Much to our surprise, in the upper right-hand corner of Figure || we see that Vequiv{r) for the Is state in Pb has 
a central attraction in addition to a dominant, mid-range repulsion. This is true even at normal potential strength 
where the changeover in sign occurs close to the r value at which Re ipi (r) has a sign change in Figure ^. Accordingly, 
we now believe that the low energy pion optical potential has an s state repulsion on the outside of the nucleus and 
an attraction within. Furthermore, this shows how, as the nonlocal optical potential is made stronger, the effective 
attractive part increases in depth and eventually leads to a confined nuclear state. 



IV. CONCLUSIONS 

We have shown again that pionic atom energies and wavefunctions can be calculated accurately in momentum 

space with a microscopic and nonlocal optical potential. Our particular potential includes effects from finite t:N 
ranges in both the first and second order terms, covariant kinematics, Pauli exclusion, nucleon recoil, off-energy shell 
dynamics, and vacuum polarization. The Coulomb force is included exactly, and in our approach the energies and 
widths are determined by searching for the poles of the T matrix for the combined Coulomb plus nuclear potential. 
With present techniques we reproduce analytic results to five significant figures. Although computations in momentum 
space use more computer time than those in coordinate space, they can handle exactly various types of nonlocalities 
and can use theoretical results with fewer approximations; we also find them more elegant. Given modern computing 
techniques and machines, widespread momentum space applications are now practical as shown by our brief survey 
of levels throughout the periodic chart. 

The focus of our calculations were the deeply bound (Is, 2s, and 2p) levels in lead which appear as hybrid Coulomb- 
nuclear states. Our calculations show these levels to be more bound and wider than reported previously : 7% 
deeper, ~ 20 — 200% wider s states, and 50 — 300% wider p and d states. The wavefunctions appear similar to 
those of Toki et al. D so we do not expect major changes in atomic formation rates. Most importantly, the states 
remain quasi-stable (nonoverlapping) and therefore are still experimentally observable. There may, however, be some 
changes needed in the experimental search for them. 

The model dependence we have found appears to arise from our inclusion of energy dependences in the optical 
potential, complex energy extrapolations, and particularly our treatment of the finite-range annihilation term. The 
use of a relativistic wave equation is important for these states, and our use of the optical potential in a relativistic 
Schrodinger equation — as opposed to the approximate Klein-Gordon equation which ignores the quadratic V^p^ term — 
may make some differences for these deeply bound states. 

We view the major uncertainty in this work to be the description of pion annihilation. We use a quasi-deuteron 
absorption model, but it is not capable of predicting precise values for the annihilation strengths Bq and Co at 
threshold, or their energy dependences below threshold. In addition, there are unknown isovector strengths to include 
for a nucleus such as lead, and they also appear to be accessible only phenomenologically and at threshold. A 
further concern than just the values of the strengths, is the simple assumption that the annihilation potential has 
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dependence; annihilation is a complicated process in which various numbers of nucleons enter, and the dependence 
is only approximately true if there were only two. 

We have also looked for pure nuclear bound states of pions at much deeper energies. We did not find any at normal 
densities and had to increase the strength of the nuclear potential at least eight-fold before any appeared. In the 
process of looking for these states we discovered that the real part of the effective local potential for a Is pion in Pb, 
while predominantly repulsive, has an inner attractive part. It is this attractive part which binds the pion within the 
nucleus as the potential strength is increased. Observing this dual character of the low energy potential experimentally 
would be fascinating, but difficult; measuring pion scattering from nuclei as a function of energy would not help much 
since the nN interaction itself is highly energy dependent. Likewise, the attractive p wave potential in Ca appears to 
have an inner part which is repulsive, and this too appears to be a new finding. 
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FIG. 1. Comparison of optical potential predictions and measurements of shifts (top) and widths (bottom) for the indicated 
shallow levels of pionic atoms. Data are from references [^- ^o|] . 

FIG. 2. Energy levels (top) and widths (bottom) of various deeply bound 7r~'^°*Pb states. The solid curves derive from the 
present p-space potential, the others from the r-space potentials of Toki et al. and Nieves et al. jl^] . 

FIG. 3. Momentum space wavefunction ipi{p) and coordinate space wavefunctions tpi{r) for n~ Pb Is state. Dashed curves 
derive from the EM potential alone, solid curves also contain the optical potential, and the dot-dashed curve shows Im ipi{r) for 
the total potential. A logarithmic scale is used for ipi{p) to show details of nuclear effect. In lower part of figure the exclusion 
of wavefunction from the nucleus of radius is ~ 7.1 fm is evident as is the node introduced into the wavefunction. 

FIG. 4. Probablity density |w;(r)|^ = \ripi{r)\^ for the Is, 2s, and 2p tt' Pb states. Dashed curves contain only the EM 
potential, solid curves EM plus optical potentials. The exclusion of wavefunction from the nucleus of radius is ~ 7.1 /m is 
evident, although the node in the real part of the wavefunction at ^ 2.5 evident in Figure ^ is not. 

FIG. 5. The diagonal momentum space potentials V{p,p) (left) and equivalent local coordinate space potentials Vequivir) 
(right) for the Is state in ■^''**Pb (top) and the 2p state in ''"Ca (bottom). 

FIG. 6. Probablity density |«;(r)p = \ripi{r)\^ of the Is 7v~ Pb wave function arising from an optical potential increased, as 
in (^3[), by 1 (solid), 8 (dot-dashed), and 15 (dashed) times its normal value. The 15-fold increase leads to a nuclear bound 
state within the nucleus of radius ~ 7.1 fm. 



TABLE I. Parameters used in modified Wood-Saxon densities {R and a in fermis). 

Rp Op Wp Ra On Wn 

'"O 2.608 0.513 -0.051 2.608 0.461 -0.051 

^^O 2.634 0.513 0.0 2.72 0.447 0.0 

*°Ca 3.669 0.584 -0.102 3.669 0.584 -0.102 

**Co 3.70 0.55 0.0 3.82 0.505 0.0 

^°^Ag 5.32 0.52 0.0 5.47 0.473 0.0 

^°^Pb 6.624 0.549 0.0 6.624 0.549 0.0 

'^°^Bi 6.609 0.545 0.0 6.88 0.5 0.0 
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TABLE II. Effect of different treatments of relativity on point Coulomb binding energies for tt~ Ca atom. Compared are 
the analytic Schrodinger equation, the numeric Schrodinger equation, the relativistic Schrodinger equation, the Klein-Gordon 
equation with reduced mass, and the Klein-Gordon equation with pion mass. All energies in KeV. 



nl 


EssiAnal) 


EsE{Num) 


Erse 




Ekge{h) 




EKGE{m^) 


Is 


-1.4809 




-1.4810 


-1.5161 




-1.5234 




-1.5287 


2s 
2p 


-0.3702 
-0.3702 




-0.37025 
-0.37024 


-0.3761 
-0.3714 




-0.3772 
-0.3717 




-0.3785 


TABLE III. Effect of different treatments of relativity on point Coulomb binding energies (in 
via Schrodinger equation, relativistic Schrodinger equation (numeric), and Klein-Gordon equation 
Za> i). 


MeV) for tt" 
(which has no 


2°®Pb atom 
s states for 


nl 






EsE 




Erse 






Ekge 


Is 






-24.97 




-42.53 






n 


2s 
2p 






-6.243 
-6.242 




-9.224 
-6.584 






n 

-6.607 


3s 
3p 
3d 






-2.775 
-2.775 
-2.774 




-3.659 
-2.921 
-2.824 






n 

-2.929 
-2.829 


4s 

4p 
4d 
4/ 






-1.562 
-1.561 
-1.562 
-1.560 




-1.930 
-1.632 
-1.592 
-1.574 






n 

-1.636 
-1.593 
-1.577 



TABLE IV. Theoretical and experimental shifts and widths in KeV. 



state 


Rcf 


ReE 


^Ipott 


^exp 


^Ipott 




Is^'^O 




25 


1 


-219.98 


-14.57 


-15.4 ± 0.1 


8.14 


7.92 ± 0.32 


2p^'^0 








-59.06 


13.4x10"^ 


14.8 ± 1.6 xlO"^ 


5.5 


6.8 ± 0.4 xlO"^ 


Is 




i 




-215.76 


-18.87 


-19.9 ± 0.1 


7.40 


6.33 ± 0.43 


2p'*°Ca 




27 




-374.60 


1.86 


1.86 ± 0.08=^ 


1.62 


1.62 ± 0.11 


2p**Ca 




i 


1 


-374.78 


1.92 


1.58 ± 0.02 


1.84 


1.60± 0.07 


3dio»Ag 




i 


1 


-922.07 


2.27 


1.97 ± 0.03 


1.24 


1.41± 0.05 


4/208pb 




i 


1 


-1582.19 


2.28 


1.67 ± 0.02 


1.1 


1.17± 0.05 


4/ 209Bi 






1 


-1622.21 


2.53 


1.83 ± 0.06 


1.31 


1.24± 0.14 



^Average experimental values. 
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TABLE V. Theoretical energies and widths in KeV for tt — Pb. The subscripts rl, r2 and r3 denote r-space calculations 
of Toki et al. Nieves et al and Konijn et al. 



nl 


Eem 


Eipott 


Erl 


Er2 


ErS 


^Ipott 


r,i 


r,2 




Is 


-12302 


-7240 


-6959 


-6778 


-6924 


754 


632 


409 


63 


2s 


-4492 


-3043 


-2962 


-2902 


-2954 


208 


183 


140 


13 


2p 


-5957 


-5399 


-5162 


-5105 


-5138 


619 


410 


275 


154 


3s 


-2262 


-1682 


-1642 


-1613 


-1633 


88 


78 


65 


5 


3p 


-2707 


-2488 


-2418 


-2395 


-2408 


201 


151 


99 


52 


3d 


-2831 


-2916 


-2854 


-2858 




125 


91 


56 




4s 


-1352 


-1067 


-1045 


-1026 




45 


40 


35 




4p 


-1538 


-1436 


-1408 


-1394 


-1394 


90 


71 


46 


24 


Ad 


-1594 


-1639 


-1606 


-1606 




73 


52 


31 




4/ 


-1579 


-1582 


-1575 


-1582 




1.1 


1.0 


1.0 
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